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We consider the Two Higgs Doublet Model (2HDM) of type III which leads to Flavour Changing 
Neutral Currents (FCNC) at tree level in the leptonic sector. In the framework of this model we can 
have, in principle, two situations: the case (a) when both doublets acquire a vacuum expectation 
value different from zero and the case (b) when only one of them is not zero. In addition, we 
show that we can make two types of rotations for the flavor mixing matrices which generates four 
types of lagrangians, with the rotation of type I we recover the case (b) from the case (a) in the 
limit tan/3 — > oo, and with the rotation of type II we obtain the case (b) from (a) in the limit 
tan/3 — > 0. Moreover, two of the four possible lagrangians correspond to the models of types I and 
II plus Flavor Changing (FC) interactions. The analitical expressions of the partial lepton number 
violating widths F (/i — > eee) and T(fj, — > ej) are derived for the cases (a) and (b) and both types of 
rotations. In all cases these widths go asymptotically to zero in the decoupling limit for all Higgses. 
We present from our analysis upper bounds for the flavour changing transition y, — ► e, and we show 
that such bounds are sensitive to the VEV structure and the type of rotation utilized. 



I. INTRODUCTION 

Flavor Changing Neutral Currents (FCNC) are forbid- 
den at tree level in the Standard Model (SM). However, 
they could be present at one loop level as in the case 
of b -> S7 |, K° -> §, K° - K° |, t c 7 

H etc. In general, many extensions of the SM permit, 
however FCNC at tree level. The introduction of new 
representations of fcrmions different from doublets pro- 
duce them by means of the Z-coupling In addition, 
they are generated at tree level in the Yukawa sector by 
adding a second doublet to the SM Q. Such couplings 
also appear in SUSY theories without R-parity . The- 
ories with FCNC were previously considered unatractive 
because they were strongly constrained experimentally, 
especially due to the small Kl — K$ mass difference. 
Nevertheless, nowadays it is hoped to observe such phys- 
ical processes in laboratory, as a result many theories 
were proposed (see above). 

Owing to the continuous improvements in experimen- 
tal accuracies, Lepton Flavor Violation (LFV) has be- 
come a very important possible source of new physics. 
Experiments to search directly for LFV have been per- 
formed for many years, all with null results so far. Ex- 
perimental limits have resulted from searches for K f — > 
fi+e- §, *> + e~ §, K+ -> Tr+M+e" |fj, 

(i+ e+7 fujl, fi+ e+e+e" @ and fi'N -> e~N 
& 

There are several mechanisms to avoid FCNC at tree 
level. Glashow and Weinberg E3] proposed a discrete 
symmetry in the Two Higgs Doublet Model (2HDM) 
which forbids the couplings that generate such rare de- 
cays, hence they do not appear at tree level. Another 
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possibility is to consider heavy exchange of scalar or pseu- 
doscalar Higgs fields (l5) or by cancellation of large con- 
tributions with opposite sign. Another mechanism was 
proposed by Cheng and Sher arguing that a natural value 
for the FC couplings from different families should be of 
the order of the geometric average of their Yukawa cou- 
plings 0. 

Taking this natural assumption and since Yukawa cou- 
plings in the SM vary with mass, it is plausible that 
the same occurs for FC couplings. Hence it is expected 
that FCNC involving the third generation can be larger, 
while the ones involving the first generation are hoped 
to be small JOI, ]i7[ . Another clue that suggests large 
mixing between the second and third generation in the 
charged leptonic sector, is the large mixing between sec- 
ond and third generation of the neutral leptonic sector. 
This is predicted by experiments with atmospheric neu- 
trinos fl8|| . 

The increasing interest in LFV processes is due to 
the strong restrictions that experiments have imposed 
on them. This consequently determines small regions of 
parameters for new physics of any theory beyond the SM. 
Some specific decays have been widely studied within 
the framework of supersymmetric extensions, because in 
Supersymmetric theories the presence of FCNC induced 
by R-parity violation generates massive neutrinos and 
neutrino oscillations |1S[| . In recent papers the decays 
fi — > &y and fi — > 3e with polarized muons have been ex- 
amined in the context of supersymmetric grand unified 
theories to get bounds in the rn~ R — \Aq\ plane pp| . 

On the other hand, a muon collider could provide 
very interesting new constraints on FCNC, for example 
/i/i — > jtxr(er) mediated by Higgs exchange [^3f which 
test the mixing between the second and third genera- 
tions. Additionally, the muon collider could be a Higgs 
factory and it is well known that the Higgs sector is cru- 
cial for FCNC |H. Finally, effects on the coupling of 
muon and tau in the 2HDM framework owing to anoma- 
lous magnetic moment of the muon could be significantly 
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improved by E821 experiment at Brookhaven National 
Laboratory f23| . 

Additionally, in the quark sector bounds on LFV come 
from AF = 2 processes, rare B-decays, Z — > bb and 
the p-parameter (2l). Reference |2l[] also explored the 
implications of FCNC at tree level for e + e~(/i + /i~) — > 

<c + tc, t -> C7(Z, gr), D° - U° and B° s -lf s . Moreover, 
there are other important processes involving FCNC. For 
instance, the decay B~{D~) — > A"~/j + t~ which depends 
on /i — t mixing and vanishes in the SM. Hence it is very 
sensitive to new physics. Another one is B~(D~) — > 
K~[/, + e~~ whose form factors have been calculated in p5|, 



The simplest model which exhibits FCNC at tree level 
is the model with one extra Higgs doublet, known as 
the two Higgs doublet model (2HDM). There are several 
kinds of such models. In the model type I, one Higgs 
Doublet provides masses to the up and down quarks, si- 
multaneously. In the model type II, one Higgs doublet 
gives masses to the up quarks and the other one to the 
down quarks. These former two models have the discrete 
symmetry mentioned above to avoid FCNC at tree level 
E3. However, the discrete symmetry is not necessary 
in whose case both doublets generate the masses of the 
quarks of up-type and down-type, simultaneously. In the 
literature, the latter is known as the model type III p5| . 
It has been used to look for physics beyond the SM and 
specifically for FCNC at tree level f2l| , In general, 
both doublets could acquire a vacuum expectation value 
(VEV), but we can absorb one of them redefining the 
Higgs fields properly. Nevertheless, we shall show that a 
substantial difference arises from the case in which both 
doublets get the VEV, and therefore we will study the 
model type III considering two cases. In the first case, 
the two Higgs doublets acquire VEV (case (a)). In the 
second one, only one Higgs doublet acquire VEV (case 
(b)). In the latter case the free parameter tan/3 is re- 
moved from the theory making the analysis simpler. 

In section II, we describe the model and define the no- 
tation we shall use throughout the document. In section 
III, we show that we can make two kinds of rotations for 
the flavor mixing matrices which generates four types of 
lagrangians, and that in the framework of the first rota- 
tion we arrive to the case (b) from the case (a) in the limit 
tan/3 — > oo, while with the second rotation we obtain (b) 
from (a) in the limit tan/3 — > 0. Furthermore, we find 
that two of the four possible lagrangians correspond to 
the models of types I and II plus Flavor Changing (FC) 
interactions. 



II. THE MODEL 

The 2HDM type III is an extension of the SM plus a 
new Higgs doublet and three new Yukawa couplings in 
the quark and leptonic sectors. The mass terms for the 
up-type or down-type sectors depends on two matrices or 
two Yukawa couplings. The rotation of the quarks and 
leptons allows us to diagonalize one of the matrices but 
not both simultaneously, so one of the Yukawa couplings 
remains non-diagonal, generating the FCNC at tree level. 

The Yukawa's Lagrangian is as follow 



ij ^CiL 

C/.07T 

ij 



o 



-4 fl i^E% (i) 



where $1,2 are the Higgs doublets, i/y and are 
non-diagonal 3x3 non-dimensional matrices and i, j 
are family indices. D refers to the three down quarks 
D = (d,s, b) T , U refers to the three up quarks U = 
(u, c, t) and E to the three charged leptons. The super- 
script indicates that the fields are not mass eigenstates 
yet. In the so-called model type I, the discrete symmetry 
forbids the terms proportional to meanwhile in the 
model type II the same symmetry forbids terms propor- 

,. i , e (7,0 D,0 E,0 

tional to £ y ' , r^- . 

In this kind of model (type III), we consider two cases. 
In the case (a) we assume the VEV as 





V1/V2 



(*2> c 





V2/V2 



(2) 



and we take the complex phase of vi equal to zero since 
we are not interested in CP violation. The mass eigen- 
states of the scalar fields are given by fl28| 



"I 



u z \ _ 
A 



h° 



cos (3 sin /3 

— sin (3 cos (3 

sin/3 \ / y/2I 

sin/3 cos/3 J y \f2Im<^\ 

cos a sin a \ / \f2Re<\>\ — v\ 

-sin a cos a J y ^/2Re<f% - v 2 



cos 
— sin 



(3) 



where tan/3 = v^/vi and a is the mixing angle of the 
CP-even neutral Higgs sector. Gz(w) are the would-be 
Goldstone bosons for Z (W), respectively. And A is 
the CP-odd neutral Higgs. i? ± are the charged physical 
Higgses. 

The case (b) corresponds to the case in which the VEV 
are taken as 



In section IV we get bounds on LFV in the 2HDM type 
III based on the decays /1 — > ej and /i — * eee. Such de- 
cays are examined in the framework of both cases (a) and 
(b) according to the classification made above, and with 
both types of rotations. We find that such constraints 
depend on whether we use cases (a) or (b) and on what 
kind of rotation is utilized. 



<*1>0 = 





vi/V2 



($2 



The mass eigenstates scalar fields in this case are 



r ± _ 

G% = 



if 



H ± = 



\/2Im</>? , A = V2Im$, 



(4) 



(5) 
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and the neutral CP-even fields are the same as in the 
former model just replacing w 2 = 0. A very important 
difference between both models is that Gz{w) is a linear 
combination of components of $1 and $2 in the model 
(a), meanwhile in the model (b) is a component of the 
doublet $1. 



III. GENERATION OF MODELS TYPE I AND 
II FROM TYPE III 

To convert the lagrangian ([!]) into mass eigenstates we 
make the unitary transformations 



D L . R = (V L , R )Dl R 
Ul,r = (T L , R )Ul R 



(6) 
(7) 



work of case (a) 



M <Mag 



AI, 



diag 



Vr 



V2 1 Vr 
Vr V2 



(8) 
(9) 



From (||), (^) we can solve for ^ D ^°^ U '° obtaining 



£D,0 = ^_ V £ M ga g y R _^ r] n,0 (1Q) 

(11) 



ffi = ^ T £ M diag TR __ ii 



V 2 



l'2 

V2 



from which we obtain the mass matrices. In the frame- 



Let us call the eqs (|l0|), (^Tj), rotations of type I, re- 
placing them into (l|) the expanded Lagrangian for up 
and down sectors are 



£ (al) = _9^P_ T j M diag KpLDH+ _ _E_U M^ 9 KP L DG + 

+ ^-^ UM?! ag U (smaH + cosa/i°) 

V2M W 8iaf3 u v ; 

--^-UM^ ag l5 UG° - l -^2lluM^ a9 l5 UA 

+ -^—Ur l u KP L DH + - — ^ UrfU [sin (a - j3) H° + cos (a - (3) h°] 

smp V2sin/3 

H — =- Ur] U j 5 UA a + h.c. + leptonic sector. (12) 

V2 sin/3 



£{ ai) = i^^U KM dia gpRDH++ 9_JJ KM di ag p RDG ^ 

{ - > Mw Mw 

j 9 7=) j r diag 



V2M w sinP 



DMp 9 D (smaH + cosa/i°) 



19 -DM% ag l5 DG° + l9 ™ t/3 DM% a9 l5 DA° 



V2M W D V2M) 



w 

1 — --Dr, t^tj+ 1 T=T„-D n [". / . q\ zjO , „„„/„. on .01 



-UKr t L 'P R DH + - DifD sin (a - /3) H° + cos (a - (3) h° 



sin/3 ' ' V2sin/3 

Dr] D ^ 5 DA + h.c. + leptonic sector. (13) 

V2sin/3 



where K is the CKM matrix. The superindex (a, /) refers 
to the case (a) and rotation type I. 

It is easy to check that if we add (|l^) and ([l3|) we ob- 
tain a lagrangian consisting of the one in the 2HDM type 
I p8|| , plus the FC interactions. Therefore, we obtain the 
lagrangian of type I from eqs ( |l2| ) and (|l3|) by setting 
i] D = rf = 0. In addition, it is observed that the case 
(b) in both up and down sectors can be calculated just 



taking the limit tan/3 — * oo. 

On the other hand, from (^), (^ we can also solve for 
yjD, o^u.o ms t ea( j f 5 to get 

AO = V?yt M *«fl Vjl _ wa f A0 (M) 

U,0 = yl T l M diag TR _^U,0 ^ 
V\ Vi 
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which we call rotations of type II, replacing them into (|l]) the expanded lagrangian for up and down sectors become 
I 



£ 



(a.,11) 
Y(u) 



■ta.nl3UM?; a9 KP L DH + - -?—UM? r ia9 KP L DG + 



M w —>---u b Mw u 

JJ M dia gjJ (< cos aH gin ah 



V2M W cos (3 
l l^lVM^ lb UA a 



V2M 



19 UM diag l5 UG 



V2M W 
1 

V2 cos (3 
-h.c. + leptonic sector 



cos (3 



U£ U KP L DH- 



U£ U U [sin {a - 0) H° + cos (a - 0) ft ] 



v^cos/? 



Ue l5 UA° 



(16) 



- £ 



(a,II) 
Y(d) 



<?tan [3 — 



UKM^PrDH- 



w 



9 UKM d ™ 9 P R DG + 



w 



_^^JL DM d ™ 9 D { cos aH° 

V2M W cos (3 y 



sin ah' 



l 9 n njdiag 



DM% a9 Tj DG° 



w 



i9tan0 -DMt ia3 l5 DA° 



1 



1 



-UK£ PrDH~* 



D !J " ' cos/3' 
D£ D D [sin (a - [3) H° + cos (a - /?) ft ] 



v^cos/? 
+ft.c. + leptonic sector 



v^cos/? 



(17) 



r 



In this situation the case (b) is obtained in the limit 
tan j3 — ► 0, for up and down sectors. Moreover, if we add 
the lagrangians ( |l2| ) and (JT^) we find the lagrangian of 
the 2HDM type II [§8) plus the FC interactions. Sim- 
ilarly like before, lagrangian type II is obtained setting 
(, D = rj u = 0. Therefore, lagrangian type II is generated 
by making a rotation of type I in the up sector and a rota- 
tion of type II in the down sector, it is valid since £ and 
£ D are independent each other and same to i] U ' D . In ad- 
dition, we can build two additional lagrangians by adding 



£ 



(o,J7) 



, . , , -r £y(<i) and ^Y(u) + £ Y(d\- S° f° ur models are 
generated from the case (a). On the other hand, interac- 
tions involving Goldstonc bosons are the same in all the 
models in the R-gauge, while in the unitary gauge they 
vanish |§. 

Finally, we can realize that in both models (a) and (b) 



J 



with both types of rotations FCNC processes vanishes 
when all Higgses are decoupled, we shall prove it by using 
the rare processes /i — > eee and /i — > ey. 



IV. LFV PROCESSES 

In the present work, we study the processes /i — > ej 
and n — > eee in the 2HDM type III. The decay width 
of /i — * e7 in both models (a) and (b) comes from one 
loop corrections, where we have used a muon running in 
the loop. The first interaction vertex is proportional to 
the muon mass and the final vertex is proportional to 
the flavor changing transition ji — > e. The decay widths 
in the two types of rotations are given by 



r^') (/x -> e 7 ) 
T (a ^ (/x - e 7 ) 



V2sin 4 /3 
\/2cos 4 /3 



| sin asm (a — /3) Fi(mjjo) + cos a cos (a — /3) F^m^o) — cos /?i<2 (w^o ) | " 



|— cos a sin (a — 13) Fi(m H a) + sin a cos (a — (3) Fi(m h o) — tan (3F 2 (m^o ) | (18) 



where 



log[:r 2 /m 2 ] 

4:TT 2 X 2 



F 2 (x) 



log [x 2 /m 2 ] 
8tt 2 x 2 



(19) 
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FIG. 2: The figure 2 illustrate the differences between the 
models (al) and (all) respect to the parameter tan (3. We have 
decoupled the higgses masses m H o = 300 GeV and m^o and 
taken a — tt/16 and m h o = 300 GeV. The curve that increases 
with f3 corresponds to the model (al). 

The decay widths for the process \x — > eee in the two 
cases read 



r(°.") (n -> eee) 



2G F mlmlN^ 
^1024^ cos 4 /3 
sin a cos (a — (3) 



cos a sin (a — (3) 



sin /3 



(20) 



And the corresponding expresions for the case (b) are 
obtained taking the appropiate limits. These FC pro- 
cesses vanish when all Higgses are decoupled. 

Now, by using the experimental upper bounds for LFV 
processes O, Q 



r( M - 



e 7 ) < 3.59 x 10~ M GeV, 



eee) < 3.0 x 10" 



GeV, 



(21) 



We see that the upper bounds imposed by \i — > e 7 are 
much more restrictive. 

We use a muon running in the loop for the calculation 
of \i — > e 7 instead of a tau as customary. This would 
be reasonable provided some conditions. If we take the 
quotient r {a > T) / lS >^ where T^'^ represents the width 
of fi — > e 7 with a muon in the loop for the case (a), 
and similarly for r( a ' T ' , and we set mjyo = 300 GeV, 
a = 7r/16 and rriA is decoupled, we can plot the quotient 



N, 



ft e 



N„ T N Te 



(22) 



by supposing that r^'^ w Y^ a ' T \ i. e., they are of the 
same order. Here N^ e denotes the FC coupling in a 
generic way. We can notice from figure 1 that the values 
obtained for the fraction cover a wide range and therefore 
this assumption is reasonable. 

We turn now to derive constraints for arbitrary values 
of the Higgs sector. Let us consider the process fi — > e 7 in 
both cases for different values of the Higgs masses and 
mixing angles. In the figure 2 we take m^o and m^a 
going to infinity. We plot N^ e vs (3 , for a — tt/16 and 
myfi = 300 GeV for the models al(all) respectively. We 
can observe that the behaviour of the models are quite 
different in a long range of tan/3. Additionally, near to 
the critical points of tan [3 the models take complemen- 
tary values. 

The 3D plots (A^ e , rrih, tua) are shown in the figure 
3 for m H = 500 GeV, a = tt/16 and tan/3 = 1. They 
represent the models (al) and (all), similar to the figure 
2. Once again, we realize that the behaviour of both 
models is quite different. 

The figure 4 corresponds to the models (all) and (bll) 
in which to#o = 300 GeV and a = tt/16. For the model 
(all) we use tan/3 = 1. These graphics illustrate that 
the cases (a) and (b) are substantially different. 



r ( Q > / ) ( M _> eee) 



2G F mlmlN^ 
^1024^3 sin 4 (3 
cos a cos (a — (3) 



sin a sin (a — (3) 



rn 

cos/3 



2 

H" 
2 



m 



h° 



in 
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V. CONCLUSIONS 

In the present work we examine a 2HDM type III which 
produces FCNC at tree level in the leptonic sector. We 
classified the model type III according to the VEV taken 



6 



ITiA 




FIG. 2 
for th 
l,the 1 




FIG. 4: Figure 4 shows the differences between models (a) 
and (b), it is plotted in the parameter space (iV Me , mu, ttia). 
The parameter tan/3 = 1 for the model (a), the higgs mass 
m H o = 300 GeV and a = 7r/16. 



by the Higgses and to the method used to rotate the 
mixing matrices. All that, in order to write down the 
lagrangian in the mass eigenstates. When both doublets 
acquire a VEV we talk about the case (a), while when 
only one doublet acquire a VEV we talk about the case 
(b). On the other hand, when we write £ D >°,£ U >° in 
terms of r) D ' Q ,r) U '° plus the mass matrices, it is called 
here a rotation of type I. Where £ D '°, £ u '° are the mixing 
matrices which couple to $2 and $2 respectively and 
YjDfi^Ufi are ^ ne pc ma t r i ces which couple to $1 and 

$1 respectively. Now, when we solve for r] D ' ,r) u ' in 
terms of £, D '°,t; U ' and the mass matrices we call it a 
rotation of type II. 

In addition, we observe that the 2HDM of type I plus 
FC interactions is generated by adding the lagrangian of 
type (a, I) in the up sector and the lagrangian of type 
(a, I) in the down sector, meanwhile the lagrangian of 
type II plus FC interactions is generated by adding the la- 
grangian of type (a,I) in the up sector and the lagrangian 
of type (a,II) in the down sector. Other two combinations 
are possible i.e. £$ff + and + £$ff . 

Moreover, if we began with a lagrangian of type (a, I) 
we would obtain the lagrangian (b,I) taking the limit 
tan/3 — > 00, while if we started with a lagrangian of 
type (a, II) we would obtain the lagrangian (b,II) in the 
limit tan/3 — ► 0. 

To illustrate the importance of this classification wc 
show graphics to find bounds on the FC coupling N^ e 
coming from the process /i — > ej and we realize that such 
bounds are sensitive to the type of rotation and also to 
the structure of the VEV. We also calculate the process 
/j, — > 3e for both kind of rotations but the constraints 
obtained were less restrictive than the ones obtained with 
the process [i — ► cy. 

Finally, to evaluate such bounds we have used a muon 
running in the loop for the process fi — > instead 
of a tau as usual. Consequently, we plot the quotient 
Nfte/ '{N^ T N eT ) in terms of m h a and (3 , getting a wide 
range of allowed values for that quotient, showing that 
this assumption is reasonable. 

We acknowledge to M. Nowakowski for his suggestions 
and for the careful reading of the manuscript. This work 
was supported by COLCIENCIAS. 



[1] M. Ciuchini, et. al, Phys. Lett. B 316, 127 (1993); Nucl. 

Phys. B 421, 41 (1994); S. Bertolini, F. Borzumati, A. 

Masiero and G. Ridolfi, Nucl. Phys. B 353 591 (1991). 
[2] H. Stern and M. K. Gaillard, Ann. Phys. 76, 580 (1973); 

C. S. Kim, J. L. Rosner and C. P Yuan, Phys. Rev. 42, 

96 (1990). 

[3] T. Imami and C. S. Lim, Prog. The. Phys. 65, 297 (1981). 

[4] J. L. Diaz-Cruz, et. al., Phys. Rev. D 41, 891 (1990); 
G. Eilam, J. Hewett and A. Soni, Phys. Rev D 44, 1473 
(1991); G. Couture, C. Hamzanoi and H. Konig, Phys. 



Rev. D 49, 293 (1995). 
[5] J. L. Hewett and T. Rizzo, Phys. Rep. 183, 193 (1989); 
G. Baremboim, et. al, Phys. Lett. B 422, 277 (1998); 
V. Barger, M. Berger and R. Phillips, Phys. Rev. D 52, 
1663 (1995); R. Martinez, J.-Alexis Rodriguez and M. 
Vargas, Phys. Rev. D 60, 077504 (1999); F. del Aguila, 
J. Aguilar Saavedra and R. Miquel, Phys. Rev. Lett. 82, 
1628 (1999). 

[6] J. Liu and L. Wolfenstein, Nucl. Phys. B 289, 1 (1987). 
[7] M. Nowakowski and A. Pilaftsis, Nucl. Phys. B 461, 19 



7 



[8] 
[9] 
[10] 

[11] 

[12] 
[13] 
[14] 

[15] 
[16] 
[17] 



[18] 
[19] 

[20] 



(1996) ; A. Joshipura and M. Nowakowski, Phys. Rev. D 

51, 5271 (1995); G. Ross and J. W. F. Valle, Phys. Lett. [21] 

B. 151, 375 (1985). 

K. Arisaka, et. al., Phys. Rev. Lett. 70, 1049 (1993). 

K. Arisaka, et. al., Phys. Lett. B. 432, 230 (1993). [22] 

A. M. Lee, et. al., Phys. Rev. Lett. 64, 165 (1990). 

R. D. Bolton, et. al., Phys. Rev. D. 38, 2077 (1983); M. [23] 

L. Brooks, et. al., Phys. Rev. Lett. 83, 1521 (1999). [24] 

U. Bellgardt et. al., Nucl. Phys. B 299, 1 (1998). 

C. Dohmen et. al. Phys. Lett. B. 317, 631 (1993). 

S. Glashow and S. Weinberg, Phys. Rev. D 15, 1958 [25] 
(1977). 

Marc Sher and Yao Yuan, Phys. Rev. D 44, 1461 (1991) [26] 
T.P. Cheng and M. Sher, Phys. Rev. D 35, 3490 (1987) [27] 

D. Atwood, L. Reina and A. Soni, Phys. Rev. D 55, 3156 [28] 

(1997) ; G: Cvetic, S. S: Hwang and C. S. Kim., Phys. 
Rev. D 58, 116003 (1998). 

Y. Fukuda, et. al., Phys. Rev. Lett. 81, 1562 (1998) [29] 
A. Kaustubh, M. Graessner Phys. Rev. D 61, 075008 
(2000). [30] 
Y. Okada, K. Okumura, Y. Shimizu, Phys. Rev. D. 61, 



094001 (2000) 

D. Atwood, L. Reina and A. Soni, Phys. Rev. D. 53, 1199 

(1996) ; Phys. Rev. D. 54, 3296 (1996); Phys. Rev. Lett. 
75, 3800 (1993). 

G. Lopez, R. Martinez and G. Munoz, Phys. Rev. D 58, 

033003 (1998] 

Marc Sher, ^ep-ph/0006159| v3 

Workshop on Physics at the First Muon Collider and at 
the Front End of the Muon Collider, ed. S. Geer and R. 
Raja (AIP Publishing, Batavia 111 1997) 
W.S. Hou, Phys. Lett B 296, 179 (1992); D. Cahng, W. 
S. Hou and W. Y. Keung, Phys. Rev. D 48, 217 (1993). 
M. Sher and Y. Yuan, Phys. Rev. D. 44, 1461 (1991). 
S. Nie and M. Sher, Phys. Rev. D. 58, 097701 (1998). 
For a review see J. Gunion, H. Haber, G. Kane and S. 
Dawson, The Higgs Hunter's Guide, (Addison- Weslwy, 
New York, 1990) 

D. Atwood, L. Reina and A. Soni, Phys. Rev. D 55, 3156 

(1997) ; Phys. Rev. Lett. 75, 3800 (1995). 
Shuquan Nie and Marc Sher, [hep-ph/9805378; 



